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Introduction 

A part of aviation accident mitigation is a crashworthy airframe structure, and an important measure of merit tor a 
crashworthy structure is the amount of kinetic energy that can he absorbed in the crush ot the structure. Prediction ot 
the energy absorbed from finite element analyses requires modeling the progressive failure sequence. Progiessive 
failure modes may include material degradation, fracture and crack growth, and buckling and collapse. The design of 
crashworthy airframe components will benefit from progressive failure analyses that have been validated by tests. 

The subject of this research is the development of a progressive failure analysis lor a textile composite, 
circumferential fuselage frame subjected to a quasi-static, crash-type load. The test data tor the frame ate reported in 
Ref. I, and these data are used to develop and to validate methods for the progressive failure response. 


Circumferential Frame Specimens and Tests 

The structural components tested in Ref. 1 were the size of a typical circumferential fuselage frame for a wide body 
commercial transport aircraft with a nominal inside radius of 1 18 inches and a radial depth ot 4.8 inches. The tiame 
segments are 48° circular arcs with an asymmetrical cross section in the shape ot the letter J. A detailed illustration 
of a frame denoted as B along with its dimensions is shown in Fig. I. 
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Figure 1 J-section fuselage frame B 
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The frames were fabricated from 2X2 2-D triaxial braided textile composite preforms coupled with the resin 
transfer molding process using 3M PR500 epoxy resin. 1 he architecture ot the pretoims is |0 /±64 39.7 '/< 

axial, and the yarns are made of AS4 graphite fibers. Material property data are listed in Table 1. 
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Table 1. Tri-axial braid properties for V f = 55*26% 


Properties 

TEXC'AD 2 

Tension 

Tests 3 

Axial modulus Ej j, psi 

7.06x1 0 h 

7.09x10° 

Transverse modulus £22^ psi 

6.59x10*’ 

N/A 

Thickness modulus E 33 , psi 

1.53x10° 

N/A 

Poisson's ratio ( v , 2 ) 

0.231 

0.26 

Poisson's ratio ( v 13 ) 

0.216 

N/A 

Poisson’s ratio ( V 2 3 ) 

0.298 

N/A 

In-plane shear modulus psi 

1.91x10° 

N/A 

Transverse shear modulus G 1 3, psi 

0.601 x! 0 ° 

N/A 

Transverse shear modulus G?3, psi 

0.645x10° 

N/A 

Tensile failure strain (ej T , ) 

14,071 

10,588 

Compression failure strain { C, , ) 

10,108 

N/A 

Axial tensile strength X 7 -, psi 

91,370 

76,880 

Axial compressive strength X c , psi 

71,000 

N/A 

Transverse tensile strength Y r , psi 

73,140 

N/A 

Transverse compressive strength Y c , psi 

56,890 

N/A 

In-plane shear strength S , psi 

30,460 

N/A 


A sketch of the frame mounted in a universal testing machine is shown Figure 2. The structural testing consisted 
of mounting frame B vertically, convex side up. in a universal testing machine and subjecting it to a slowly applied, 
radially inward displacement at its apex. In addition to in-plane bending and circumferential compression, the 
asymmetric 7-section frame exhibits out-of-plane bending and torsion. The largest circumferential strain magnitude 
was compressive and occurred in the outer flange at the apex. A compressive hoop strain of approximately one-hall 
this maximum value occurred in the inner flange near the supported ends of the frame. The largest circumferential 
tensile strains are much less than the largest compressive strain magnitudes. The maximum compressive strain 
recorded at the first major failure event is about 4000 pe which is substantially less than the material ultimate 

compression strain of 10,108 pe predicted by TEXCAD 2 computer program for this material system. 
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Figure 2 Sketch of the frame test apparatus 


Frame B Progressive Failure Analysis 

Finite element models of frame B were developed in the ABAQUS/Standard software package 4 using the dimensions 
listed in Figure 1 , the material property data listed in Table 1 , and using the finite-membrane-strain shell elements 
available in ABAQUS. Clamped boundary conditions were prescribed in the analyses to simulated the frame mounted 
in the end blocks. To simulate the contact of the platen with the frame, we used a contact algorithm available in 
ABAQUS. A rigid plane parallel to the platen was displaced in the downward direction into the outer flange of the 

frame. 

Intralaminar Failure and Degradation 

A FORTRAN subroutine was written to be used in conjunction with ABAQUS to predict intralaminar failure initia- 
tion followed by a degradation in the moduli. The algorithm follows one used in Ret. 5 tor intralamina progressive 
failure of a unidirectional composite material. The in-plane strengths for the triaxial braid determined from 1 EXCAD 
are listed in Table I. With this limited material strength data, we assumed a maximum in-plane stress criterion to pre- 
dict the initiation of failure in an element. The FORTRAN subroutine used is called from the ABAQUS input deck 
under the *USER DEFINED FIELD command. The state of stress in an element from the equilibrium state at the pre- 
vious load is passed to the subroutine prior to the next load increment. The in-plane stresses are compaied to their 
corresponding strengths. If the stress magnitudes equal or exceed their corresponding strengths, then the value ol a 
field variable is changed from zero to one. Field variables are a means to indicate whether failure has occurred or has 
not occurred; a value of zero means no failure and a value of one means failure. The values ol the field vai iables ate 
returned to the input deck where they are used to determine which lines of material properties are used to calculate the 
stress and strains for the next load step. 

The maximum in-plane stress criterion is non-interactive and is based on the values of the axial normal stress 
S, , , the transverse normal stress S 22 . and the in-plane shear stress S 12 . Consequently three field variables aie 
required: Field variable FV 1 is associated with the axial normal stress limits, FV 2 is associated with the transverse 
normal stress limits, and Held variable FV 3 is associated with the shear stress limits. The initial values of the three 

Held variables are zero, representing an intact layer. If S n >91. 37ksi or S M < -71 ksi , then FV 1 = 1 and the axial 
modulus (E n ) is multiplied by a factor of 1/I0 (> and the transverse modulus (E 22 ) and the shear moduli (G , 2 , G, h 
C 2i ) are multiplied by a factor of 1/1000. If S 22 > 73. 14ksi or S 22 < -56.89ksi , then FV 2=1, and the same 
scheme as used before is followed but the transverse modulus (E 22 ) is multiplied by a factor of 1/10 6 and the axial 


3 of 8 


Grant NAG -I -Oil 23 




modulus (£/,) is multiplied by a factor of 1/1000. If |S I2 | > 30.46ksi , then FV 3 = 1 and the in-plane moduli (E„, 
) and the out of plane shear moduli (G /} , G 2 ,i) are multiplied by a factor of 1/1000. the in-plane shear modulus 
(G /2 ) is multiplied by a factor of 1/I0 6 . For combinations of failure modes, the previous reduction factors are 

combined together decreasing each property by the largest factor unless the factors are equal and greatei than 1/10 . 
In the latter case, there is an additional reduction of 10. For example, it the values of field variables aie (1,0,1), which 
indicates a combination of axial and shear failure, then E n and C«| 2 are decreased by 1/10 6 and H 22 . G l3 , G 23 are 
multiplied by 1/10 4 because the axial failure factor is I/10 3 and the shear failure factor is 1/I0 3 . See Table 2 for the 
property degradation scheme. 


Table 2. Reduction Factors For Intralaminar Material Degradation 


ABAQUS field variables 

Reduction Factor for Each Material Property 

FV 1 

FV2 

FV3 

Eli 

E 22 

e 33 

V j 0 

v 13 

v 23 

G,2 

G,3 

G 23 

0 

0 

0 

I 1 1 1 1 1 1 1 1 

1 

0 

0 

I/10 6 

I/10 3 

1 

0 

0 

0 

1/10 3 

1/10 3 

1/I0 3 

0 

l 

0 

1/I0 3 

l/IO 6 

1 

0 

0 

0 

I/10 3 

1/10 3 

1/I0 3 

0 

0 

1 

1/1 o 3 

I/10 3 

1 

0 

0 

0 

l/IO 6 

1/1 o 3 

1/10 3 

1 

1 

0 

l/IO 6 

l/IO 6 

1 

0 

0 

0 

l/IO 4 

l/IO 4 

l/IO 4 

1 

0 

1 

l/IO 6 

l/IO 4 

1 

0 

0 

0 

1/10 6 

1/10 4 

l/IO 4 

0 

1 

1 

l/IO 4 

1/10* 

1 

0 

0 

0 

l/IO 6 

1/10 4 

l/IO 4 

1 

1 

1 

l/IO 6 

1/I0 6 

1/I0 3 

0 

0 

0 

1/10 6 

1/10 5 

l/IO 4 


Delamination Initiation and Progression 

Interface elements are used to model the delamination process between adjacent layers of the quasi-laminar textile 
composite. The interface element used ter model interlaminar crack initiation and growth in tape laminates developed 
by Goyal, Johnson, Davila, and Jaunky 6 was implemented in the finite element analysis (FEA) of the frame. Interface 
elements were located in the models at the mid surface of the outer flange, where delamination was observed to ini- 
tiate in the tests. The displacements of the sublayers above and below the mid plane where the interface elements aie 
located determine interlaminar tractions through a phenomenological constitutive law. Each interface element con- 
sists of two surfaces and eight nodes, an upper and lower surface with four nodes each. The tour nodes on the uppei 
surface of the interface element connect to the four nodes on the lower surface of the element in the upper sublayer of 
the outer flange. The four nodes on the lower surface of the interface element connect to the four nodes on the upper 
surface of the element in the lower sublayer of the flange. In the undeformed state the two surfaces of the interface 
element coincide. The constitutive law that resists the separation of the upper and lower surfaces of the interface ele- 
ment can be conceptualized as representing a continuous distribution of nonlinear springs that do not interact with 
one another; i.e. similar to Winkler foundation. A graphical representation of the interface element is shown in 
Figure 3. 

The constitutive law relates the tractions to the displacement jumps across the interface. An exponential form of 
the this relationship is used in Ref. 6, which is illustrated for Mode I loading, the opening mode, in Figure 4. In this 
figure T is the traction, A is the displacement jump normal to the interface, T‘ is the specified tensile, or peel, 
strength, and G, is the specified critical strain energy release rate in Mode I. Tractions first increase with increasing 
A . reach a maximum at A = A, , then decrease as A > A, . The exponential form of the law and specification of T' 
and G, determine the critical separation length A = A,. . The critical separation length is a characteristic length scale 
for the initiation of two free surfaces of the incipient crack. The softening portion of the constitutive law repiesents 
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Figure 3 Upper (S + ) and lower (S') surfaces of the interface element. 


the presence of the process zone ahead of the crack tip, which is a craze zone in polymers or a plastic zone in ductile 
metals. Under combined Mode I, II, and III loading, delamination initiation is predicted using multi-axial stress 
criterion and the propagation of delamination is based on an mixed-mode fracture criterion 6 . A damage parameter is 
included in the formulation such that on unloading from a state where A > A,. . the traction does not follow the loading 
curve. That is, unloading is associated with a loss of work. 
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Figure 4 Representative traction-stretching curve for the springs 


Because there is little test data for debonding strengths and fracture energies of 2D triaxial braided textile 
materials, and testing for this data is beyond the scope of the project, values used in the analyses are based on those 
used for unidirectional graphite epoxy composite laminates in Ref. 6. Values used tor the three interfacial strengths 

T< , / = 1 , 2, 3 , and the three critical strain energy release rates G lr , G //r , G llh . are listed in Table 3. 


Table 3. Interfacial Strengths and Fracture Energies from Ref. 6 


T\ , T‘ 2 

T\ 

G/r 

G//, > G ///r 

1 1.6 ksi 

8.7 ksi 

2.01 Ibs/in 

8.28 lbs/in 
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Results 

Refinements in the forty-eight degree frame model were required to capture the delamination process in the outer 
flange. To keep the model of reasonable size for the computational resources available, the model size had to be 
decreased for the smaller elements required to model decohesion associated with delamination. Using symmetry, only 
half of the frame was modeled. Boundary conditions modeling symmetry were used at the midspan location. With the 
physical model size decreased, the element size was reduced so that the smallest element in the outer flange was on 
the order of 0.06 inches (~1.6 mm). The model consisted of S4R shell elements. R3D8 rigid elements for the platen in 
the contact analysis, and the user-defined interface elements. The final model contained 16.770 elements and 15,08 1 
nodes. However, ABAQUS then created 2,268 additional elements and 6,195 additional nodes of it’s own to run the 

contact analysis. 

Interface elements are only included in the midplane of the outer flange over a ten degree arc beginning at the 
midspan because of modeling difficulties and size limitations. Hence, the model cannot represent crack propagation 
from the outer flange into the web at the apex of the frame. Intralaminar progressive failure prediction is included tor 
the shell elements making up the outer and inner sublayers of the outer flange surrounding the interface elements, and 
for the elements comprising the central span of the web and a section of the inner flange near the clamped end. This 
analysis captured the load-shortening response behavior fairly well with a large load drop followed by reloading as is 
shown in the load-deflection plot of Figure 5. The slope of the reloading portion of the response from FEA is within 
10% of the lest data at reloading. The energy absorbed predicted by the FEA is 11% lower than measuied by the test, 
where the energy absorbed is the area under the load-deflection plot. 

Along with the good correlation of the load-deflection behavior, the location of element failures correlate well 
with the locations of failure observed in the test. A comparison of the tailed frame and the predicted element lailures 
is shown in Figure 6. Failures observed during the test of frame B included the separation of the axial filler material in 
the outer flange at the apex of the frame, cracking of the front and back sides of the outer flange, and cracking of the 
web from the outer to the inner flange. Delamination and intralaminar element failures were predicted at the apex of 
the frame over the junction between the web and outer flange. Intralaminar failures were predicted in the back side 
outer flange and the web near the apex. The failure pattern and sequence predicted by the progressive failure model 

correlated very well with damage observed during the failure of frame B as reported in Ref. 3. 


Concluding Remark 

Including both interlaminar and intralaminar progressive failure models, geometric nonlinearity, and the contact of 
the platen with the frame in the ABAQUS finite element analysis resulted in very good correlation between the analy- 
sis and test through several major failure events of Irame B. 
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Displacement of Mid Point (in) 

Figure 5 Load displacement responses from the analysis including both the intra- and inter- 
laminar progressive failure models and the test. The domain of the intralaminar failure 
extended to web and inner flange. 
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Figure 6 Comparison of the element failures predicted by the 
progressive failure model and failures that occurred on frame B. 
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